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Abstract. In this paper, I compute the inhomogeneous (multi-probability) bond 
critical surfaces for the (4,6, 12) and (3 4 ,6) lattices using the linearity approximation 
described in (Scullard and Ziff, J. Stat. Mech. P03021), implemented as a 
branching process of lattices. I find the estimates for the bond percolation thresholds, 
p c (4,6, 12) = 0.69377849... and p c (3 4 ,6) = 0.43437077..., compared with Parviainen's 
numerical results of p c « 0.69373383 and p c ps 0.43430621 . These deviations are of 
the order 10 -5 , as is standard for this method, although they are outside Parviainen's 
typical standard error of 10~ 7 . Deriving thresholds in this way for a given lattice 
leads to a polynomial with integer coefficients, the root in [0, 1] of which gives the 
estimate for the bond threshold. I show how the method can be refined, leading to a 
sequence of higher order polynomials making predictions that likely converge to the 
exact answer. Finally, I discuss how this fact hints that for certain graphs, such as the 
kagome lattice, the exact bond threshold may not be the root of any polynomial with 
integer coefficients. 
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1. Introduction 

Despite its deceptively simple definition, percolation provides a wealth of interesting 
problems that have kept physicists and mathematicians busy for over fifty years pQ. 
In bond percolation, which is the focus of this paper, we declare each bond of a 
lattice to be open with probability p and closed with probability 1 — p, resulting in 
a random distribution of connected clusters. In the limit of an infinite lattice, there 
is a sharp critical probability, p c , at which an infinite cluster first appears. The 
determination of p c is an unsolved problem except for a few 2-dimensional cases. There 
are also many other quantities of interest, and recently, arguments from conformal 
field theory [2] and the discovery of Schramm-Loewner evolution [3J H] have facilitated 
some spectacular calculations in the continuum limit, in which details of the underlying 
lattice disappear ([3 EJ [7] are just a few examples). However, many new lattice- level 
results have also appeared in recent years, including exact calculations of percolation 
thresholds [HI El HDl EU CGI EG] as well as other rigorous [El |15j [16] and numerical 
[T7J [HI [191 [20] results at the critical point. This paper describes an approximation 
scheme that allows very accurate determination of bond percolation thresholds for 
arbitrary periodic lattices. It was originally discussed in [21] . and expanded in [22] 
where thresholds were estimated for all but two of the Archimedean lattices to within 
10 -5 of the numerically determined values. I complete this program here by computing 
the approximate thresholds for the (4, 6, 12) and (3 4 , 6) lattices. In addition, I show 
that the approximations can be refined, leading to a sequence of polynomials that give 
successively better predictions for the bond thresholds. Finally, I discuss how this 
refinement provides significant clues about what the eventual solutions of these problems 
might look like. 

2. Critical polynomials 

Every two-dimensional lattice for which the percolation threshold is known exactly is 
in the category of 3- uniform hyperlattices [13]. An example of such a hyperlattice is 
shown in Figure [3j The unit cell of the graph is contained between three vertices, but 
the shaded triangle can represent any connected network of sites and bonds. The critical 
points of all these problems are given by the Ziff criterion [9], 

P{A,B,C) = P{A,B,C), (1) 

where P(A, B, C) is the probability that all vertices (A, B, C) can be connected within 
a shaded triangle, and P(A, B, C) is the probability that all three are disconnected. 
Thus, the critical threshold, which marks the appearance of an infinite connection, can 
be located by comparing probabilities of events on a single unit cell. This is a very 
special property that is not shared by most lattices. The result is that in this class, 
every bond threshold, for example, is the root in [0, 1] of a polynomial of degree at 
most n with integer coefficients, where n is the number of bonds in the unit cell. That 
is, at present, all exactly known thresholds are algebraic numbers. The criterion (CQ) 
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Figure 2. The (3 4 ,6) lattice with the assignment of probabilities on the unit cell. 



provides the opportunity to derive not only critical thresholds, but also critical surfaces 
for inhomogeneous percolation (also called anisotropic percolation in some applications 
[23J). Here, each bond, i, of the unit cell has its own probability, pi, and applying (JT]) 
leads to a condition of the form 

/( Pl ,p 2 ,..,p n ) = (2) 

where none of the pi appears as a power greater than one. For example, for the square 
lattice, with probability pi on the vertical bonds and P2 on the horizontal bonds (Figure 
B|a)), we have, 

S(pi,P2) = 1 -Pi -P2 = (3) 
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and for the 3-probability honeycomb lattice (Figure IH(b)), 

H(pi, P2, Pa) = 1 - P1P2 - P2P3 - P1P3 + P1P2P3 = 0. (4) 

A more complicated example is the martini lattice |8] El], which is derived in [22] and 
shown in Figure fflc), 

M(pi,p 2 , p 3 , ri, r 2 , r 3 ) = 1 - p\p 2 r z - P2P3^i - Pmr 2 - P\P2r\r 2 

- PlP3^1^3 - P2P3^2^3 + V\V2Vzr\T2 

+ PiP2P3^i^3 + PiP2Psr 2 r 3 + Pip 2 r ± r 2 r 3 

+ pmr^r-s + p 2 P3r 1 r 2 r 3 - 2p 1 p 2 p 3 r 1 r 2 r 3 = . (5) 

Setting the probabilities equal in these equations gives the corresponding critical 
polynomial, whose root in [0, 1] determines the threshold. A question that arises is 
whether other lattices, not in the 3-uniform class, also have algebraic critical points. 
So far, the answer to this question is not known. However, it is possible to find 
approximations to a range of unsolved bond problems by effectively assuming that the 
critical threshold is the root of a polynomial of degree n, where n is the number of 
bonds in the unit cell. This method has been used with success in [21] and [22], and 
a corresponding approximation for the Potts model was used by Wu for lattices of 
the kagome type [251 EEl ETj. Here, I use this method to find the critical polynomials 
approximating the thresholds of the (4, 6, 12) and (3 4 , 6) lattices, and I show that this 
leads to a sequence of approximations that likely approach the exact answer for any 
periodic lattice. 

I will begin by defining some terminology. All the lattices considered here are 
periodic and, as we have done already, we call the smallest graph that can be copied 
and translated to form the entire lattice, the unit cell. When we consider inhomogeneous 
percolation and assign different probabilities to different bonds, we are not confined to 
remain within a single unit cell; corresponding bonds on neighbouring cells may have 
different probabilities. In this case, the percolation process defined on the lattice does 
not have the same periodicity as the lattice itself. I will term the smallest periodic unit 
of the probabilities the "base" of the process. Examples for the (4, 8 2 ) lattice are shown 
in Figure El The lattice is in Figure |5^a), Figure [5](b) depicts a base consisting of a 
single unit cell, and the base in Figure [5]^c) has 12 different probabilities and covers two 
unit cells. 

3. Linearity assumption 

The approximation used here was described fully in [22], but I will give a brief summary. 
Note that in equations ([3]) and (@J, all probabilities appear only in first order. This is 
a natural consequence of (111), and I refer to this property as "linearity". As mentioned 
above, this is a very special property, but it turns out that very good approximations to 
bond percolation thresholds can be found by assuming that it holds in general. It may 
not be very clear why such a strategy should produce good approximations. However, 
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Figure 3. A 3-uniform hyperlattice and unit cell. The shaded area can represent any 
network of sites or bonds. 



Pi 



Pi 
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Figure 4. a) the square lattice; b) the honeycomb lattice; c) the martini lattice. 
Inhomogeneous probability assignments are shown below the lattices. 

the study in [22] shows that the predicted thresholds are generally within 10~ 5 of the 
best numerical estimates, and that continual refinement is possible. 

Consider the (4, 8 2 ) lattice (Figure |^ a)), in which each bond on the unit cell has 
a different probability, and call its critical surface FE(p, r, s, t, u, v). If we remove the 
p— bond by setting p = 0, the result should be the honeycomb lattice, with the trivial 
difference that two of the bonds are doubled. Contracting the p— bond to zero length by 
setting p = 1 gives the martini- A lattice (Figure E^a)), which is another exactly known 
threshold, denoted A(pi,p 2 , r\, r 2 , r 3 ) and given by 

A(pi,p 2 ,»"i, r 2 ,r 3 ) = M(pi,p 2 , l,r 1 ,r 2 ,r 3 ). (6) 

The only way to satisfy both of these conditions and keep the function FE(p, r, s, t, u, v) 
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Figure 5. a) the (4,8 2 ) lattice; b) inhomogeneous assignment of probabilities on the 
unit cell; c) a base that consists of 12 probabilities and covers two unit cells. 

first-order in its probabilities is to set 

FE(p, r, s, t, u, v) = pA(r, s, t, u, v) + (1 — p)H(s, tv , ru) = 0. (7) 

Expanding gives, 

FE(p, r, s, t, u, v ) = 1 — (pru + stu + psv + rtv) — {rsuv + ptuv ) + 
prsuv + ptruv + pstuv + rstuv + prstu + prstv — 2prstuv. 

The prediction for the bond threshold is found by solving the polynomial equation 
FE(p,p,p,p,p,p) = 0, 

1 - 4p 3 - 2p 4 + 6p 5 - 2p 6 = 0, (8) 

with solution on [0, 1] p c = 0.676835... . Comparing with the numerical estimate 
p c ~ 0.676802 with a standard error of 6.3 x 10~ 7 , we find that although our solution 
is very close, it is ruled out numerically. The assumption that FE is first-order in 
its arguments is therefore not correct. However, in seemingly all cases, the linearity 
assumption produces good approximations to bond thresholds. 

At this point, one may wonder if the approximation is actually well-defined or if 
the derived critical surface depends on the method used to find it. It is clear that 
equation ([71) is the only way to satisfy reduction to the honeycomb and A lattices in 
the appropriate limits of p. What is not obvious is whether this choice is consistent 
with the results one expects by setting bonds other than p to or 1. For example, 
the f-bond is not equivalent to the p-bond but we may just as well use v to derive 
the threshold. If we set this bond to 0, we disconnect the graph and end up with one- 
dimensional percolation with the decorated bond shown in Figure [6](b). The probability 
of crossing this bond is given by Pd{p, r , s, t, u) = u[l — (1 — pt)(l — rs)}, and, since this 
is 1-D percolation, our critical surface should be pr> = 1 or 1 — p D = 0. Setting v = 1 
gives the lattice shown in Figure E](c). The critical surface of this lattice, which for 
the moment I will call Li (actually, it is the dual of Wierman's bow-tie lattice [28J), is 
not known exactly for this configuration of probabilities (however, see [22]), so we must 
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a) b) c) d) 

Figure 6. Lattices and probability assignments used to derive the approximate critical 
surface ([?]) ■ a) the martini- A lattice; b) one-dimensional percolation with a decorated 
bond; c) the bow-tie dual; d) the martini-B lattice. 

employ the same procedure to find Li(p,r, s,t,u). Setting s = in this lattice gives 
us the exactly- known martini-B lattice [HI E] of Figure [6](d) [HI |9], with critical surface 
B(pi, n, r 2 , r 3 ) = M(pi, 1, 1, ri, r 2 , r 3 ), while s — 1 is once again the honeycomb lattice 
with a doubled bond. Thus we have, 

Li(p,r,s,t,w) = sB(u,t,p,r) + (1 - s)H(p,r,ut), (9) 

and now 

FE(p, r, s, t, u, v) = vLx(p,r, s,t,u) + (1 - v)[l - p D (p, r, s, t, «)] = 0. (10) 

Expanding gives the same function as in (171) . The answer is therefore unchanged whether 
we choose to impose constraints on the behaviour of p or v. We would like to know 
whether this will be true in general. Given a periodic lattice with a base of n bonds, 
can we find a critical function, first order in all pi, such that the results of setting each 
of the pt alternately to and 1 are consistent with the answers required in those limits? 
To answer this question, we observe that the linearity assumption is equivalent to the 
requirement that the critical function satisfy the n— dimensional Laplace equation in the 
probabilities. Specification of the limiting functions when Pi = and Pi = 1 for all i is 
equivalent to the imposition of boundary conditions on the faces of an n— dimensional 
hypercube. Provided the boundary conditions are consistent with one another, existence 
of the function is assured by the standard existence theorem for the n— dimensional 
Laplace equation. This also proves uniqueness, although it was already clear from 
the derivation of (jTJ). Thus, we can talk about the first-order function that serves as 
the approximate critical surface for a given base, and the method we use to find it is 
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irrelevant. Critical functions were found in [22] for all but two of the Archimedean 
lattices by starting with the most general function first-order in its arguments, and 
imposing symmetry and special cases until all unknown constants were fixed. However, 
the greater the number of bonds in the base, the more cumbersome this procedure 
becomes, and it was not possible to find the thresholds of the (4, 6, 12) and (3 4 , 6) 
lattices in this way, which have 18 and 15 bonds in their unit cells. In this paper, I 
will find these two critical surfaces using the technique outlined above for the (4, 8 2 ) 
lattice; choose a bond, call L the lattice resulting from setting Pi = for some i, and 
Li the lattice resulting from setting p^ = 1. If the threshold of either of these lattices, 
say L , is unknown then repeat the process on that lattice, finding L 00 and L m . By 
continuing in this way, the number of bonds will eventually be reduced sufficiently that 
lattices with exactly known thresholds appear. Although the result can be a fairly large 
branching process of lattices, the procedure is straightforward and, aside from some 
trivial computer algebra, requires minimal computational resources. 

4. Archimedean lattices 

4.1. (4,6,12) lattice 

The (4, 6, 12) lattice is shown in Figured] and we denote its threshold by 

FST(pi,p2,P3,P4,P5,P6,ri,r2,r3,r4,r 5 ,r 6 ,si,S2,ii,^2,Wi,W2)- (11) 

The unit cell contains 18 bonds, so we have no hope of encountering a known lattice until 
we have performed many iterations. We begin the journey by setting p\ — 1, resulting 
in the graph L\ (Figure [TlTa)) with a bond doubled in series, and setting pi = 0, giving 
L (Figure [TOT a)). Their critical surfaces are given by the functions 

L (P3,P4,P5, ri, r 2 , r 3 , r 4 , r 5 , r 6 , s u s 2 , t 2 , u u x, y) (12) 

and 

£i(P2, P3, P4, Ps, Pa, r x , r 2 , r 3 , r 4 , r 5 , r 6 , s u s 2 , h, t 2 , u u u 2 ) (13) 

and thus we have 

FST = (1 - pi)L {p 3 ,p A ,p 5 ,ri,r 2 ,r 3 ,r4,r 5 ,r 6 ,s 1 ,s 2 ,t2,u 1 ,tip 2 ,u 2 p 6 ) 

+ piLi(p 2l p^ p 4 , p 5 ,p e , n, r 2 , r 3 , r 4 , r 5 , r 6 , s u s 2 , t u t 2 , m, u 2 ). (14) 

Since neither of these functions are known, we start the procedure over, beginning with 
Lq. The rest of the gory details are described in the Appendix and all the lattices arising 
in the branching process are shown in Figures \IU\ and [TT] with probability assignments 
in Figures [121 and [TBI The final critical function TFS has 1932 terms, and is included in 
a text file in the supplementary material to this submission. The polynomial resulting 
from setting all probabilities equal is 

1 - 18p 6 - 6p s + 30p 9 + 3j9 10 + 108p n - 81p 12 - 174p n 
- 246p 14 + 1090p 15 - lllOp 16 + 480p 17 - 78p 18 = (15) 
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with root in [0, 1], p c = 0.69377849... . According to the numerical work of Parviainen 
[29J, p c = 0.69373383..., putting our prediction within 4.5 x 10~ 5 of the numerical value, 
but outside his standard error of 7.2 x 10~ 7 . This level of accuracy is typical for the 
method [22]. 

Before arriving at this answer we encountered many intermediate lattices (there 
are 21 in Figures ITUl and fTTT) . and the procedure may rightly be described as tedious. 
However, the upside is that we obtain predictions for all these lattices as a by-product, 
some of which are even exact. Table [1] shows the predictions along with an indication 
of which are exact. 

4.2. (3 4 ,6) lattice 

The procedure for the (3 4 , 6) lattice (Figure |2]) is described in the Appendix. The full 
critical surface has 12795 terms and is included in the supplementary material. The 
polynomial is 

1 - 12p 3 - 36p 4 + 21p 5 + 327p 6 - 69p 7 - 2532/ + 6533p 9 - 8256p 10 
+ 6255J9 11 - 2951j9 12 + 837p 13 - 126p 14 + 7p 15 = 0. (16) 

This predicts p c = 0.43437077..., whereas Parviainen gives p c = 0.43430621..., so again 
we are fairly close, differing by 6.5 x 10~ 5 but still outside Parviainen's standard error 
of 5.0 x 10~ 7 . The predicted thresholds for the intermediate lattices are shown in Table 

El 
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Table 2. Predicted thresholds for the lattices in Figures [14] and [l5l(*) denotes exact 
threshold. Graphs with threshold 1/2 are self-dual. 
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5. Polynomial sequences 

Every critical bond threshold derived in this way is the root in [0, 1] of a polynomial 
of degree at most n, where n is the number of bonds in the base. However, as already 
mentioned, the base need not consist of only a single unit cell. Considering the base of 
the (4, 8 2 ) lattice to be two unit cells results in a 12 th order polynomial, which I will 
call the second polynomial, that makes a better prediction (see below) than that found 
by the first polynomial from the 6— bond case. It is natural, then, to assume that by 
considering larger and larger bases, one can get arbitrarily close to the exact threshold. 
For those lattices to which ([1]) applies, the first polynomial gives the exact answer. For 
those situations, all subsequent polynomials must also give the same result even if the 
full critical surface is not correct. 

One possible mathematical question is whether there is always only one root in [0, 1] 
of the polynomial when following the above procedure. I do not know of any argument 
that establishes this in general, but for all lattices I have encountered, it has always 
been true. 

5.1. Second polynomials 

The polynomials calculated above represent only the first step in a series of 
approximations. By extending the base to include a second unit cell, we find a higher- 
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Figure 7. a) the kagome lattice; b) the unit cell; c) the base using two unit cells. 

order polynomial for the critical threshold. This generally, but not always, gives a 
refinement of the estimate and I explore this idea for a few lattices. 

5.2. (4,8 2 ) lattice 

We already encountered the first polynomial for this lattice, but we can get a better 
estimate using the base shown in Figure 0(c). The full threshold can be found in the 
supplementary material, but the critical polynomial is 

1 - 4p 4 - 16/ + I2p 7 + 22p 8 + 16p 9 - 70p 10 + 48p n - 10p 12 = (17) 

with p c = .676787... . The difference between this and the numerical value is 1.5 x 10~ 5 
as opposed to 3.3 x 10~ 5 for the 6— bond estimate, so we have cut our error in half. 
Furthermore, setting corresponding probabilities on each cell equal, pi = P2 = P, 
ri = r 2 = r, etc., we get a prediction for the 6— bond case that does not reduce to 
what we found before (equation ([7])). Note also that the second polynomial is not just 
the first with a few extra terms, but is completely different. 

5.3. Kagome lattice 

The kagome lattice is shown in Figure [3(a). The threshold using the unit cell as base 
(Figure [3(b)) was derived in [10] and [22], where it was found that p c = 0.524430... 
compared with the recent numerical estimate 0.524405... [20j [27]. In considering the 
12— bond base (Figure 0(c)), I found the same prediction as in the 6— bond case, i.e. 
the first polynomial factors out of the second. Similarly, setting the corresponding 
probabilities in the neighbouring cells equal, p\ = pi = p,ri = r 2 = r, etc., the 6— bond 
critical surface factors out of the larger expression. Presumably, one needs to consider 
an even larger base to get any refinement. 
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Figure 8. a) the (3,4) lattice; b) unit cell; c) a base employing two unit cells 



54. (3 3 ,4 2 ) lattice 

Both the first and second polynomials for this lattice were derived in [22], but I will 
repeat them here because they will be useful later. The (3 3 , 4 2 ) lattice is shown in Figure 
Ufa). The first polynomial, derived using the base in Figure Mjo), is 

l-2j9-2j9 2 + 3j9 3 -p 4 = (18) 

with critical threshold p c = 0.419308... whereas Parviainen gives p c = .419642... . We 
refine the estimate by considering the base consisting of two unit cells, shown in Figure 
|8]^c). In this case, the procedure gives the 10 th -order polynomial, 

1 - Ap 2 - 12p 3 + 104p 5 - 193/ + 146p 7 - 45p 8 + 2p w = (19) 

and p c = 0.419615..., which is closer to the numerical value. Once again, these two 
polynomials do not appear to have much in common other than a similar root in [0, 1]. 
However, that is not to say that no relationship exists, and in fact a recursive method 
of generating successive polynomials might even be considered an exact solution of the 
problem. 



6. Limiting case 

It is natural, after calculating a few first and second polynomials, to wonder whether 
the process must really be carried out indefinitely, or if there are some lattices for 
which an m th polynomial actually provides the exact answer. Indeed, in the case of 
the 3— uniform lattices, the first polynomial is exact and all subsequent polynomials 
make the same prediction. Could it be that in some cases we just need to consider a 
large enough base and we will find that the exact critical surface is first-order in the 
probabilities? I will argue that for many of the lattices considered here, the answer is 
no; if the first critical surface is not exact, then no finite-sized base will give the exact 
answer. 

Consider the (3 3 ,4 2 ) base shown in Figure |9]^a), where we have extended it to 
include m unit cells, and thus we have 5m probabilities, denoted by {p}. Assume also 
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that we have found the threshold by the linearity assumption, TF( m ) ({p}) , and we want 
to know if it might be exact. To answer this, we may set Pi = rj = tj = for i > 1 to 
give the lattice shown in Figure [9](b), which consists of unit cells of the (3 3 ,4 2 ) lattice 
with some long straight paths (shown in red) connecting them. Setting the probabilities 
of all the bonds in those paths to 1, we recover the case in which the base is the unit cell 
(Figure E^b)), i.e. we are left with 5 probabilities and a critical surface that is first order 
in these. The prediction for this function must agree with the one we found previously, 
as the linear threshold is unique, but that prediction has been ruled out numerically. 
Moreover, the critical surface TF^ m \{p}) is not internally consistent as it makes two 
different predictions for the same lattice. Any base consisting of a finite number of units 
cells will suffer from this defect. We therefore come to the conclusion that although 
increasing the size of the base of the process on the (3 3 ,4 2 ) lattice makes the linearity 
approximation more accurate, it will never be exact for any finite-sized base. 

For some lattices, one may be able to construct a large but finite base that does 
not suffer from this inconsistency. However, for many problems, such as the kagome 
and (4, 8 2 ) lattices, this does not seem possible, suggesting that, for these problems, the 
threshold is not the root of a finite polynomial. However, a few issues would have to be 
resolved before we could reach that conclusion. Just because the linearity hypothesis 
fails for any finite base, we cannot conclude that the exact critical surface for a base using 
only a unit cell does not contain only a finite number of terms. Consider the kagome 
critical surface for the unit cell (Figure Efb)), K(p, r, s, t, u, v). Is it possible that there is 
a maximum order, m, of the terms in the exact threshold, so that the highest order term 
possible is p m r m s m t m u m v m ? it is certainly true that if we had a base of m unit cells, and 
found the linear threshold, it would predict that K (p, r, s, t, u, v) has maximum order m. 
However, it is not clear if the converse is true. If K(p, r, s, t, u, v) has maximum order 
m, then is it necessarily derived from a critical surface that was first-order on some base 
of m cells? If we could show that this is the case, then we would know that the exact 
critical function K(p, r, s, t, u, v) is necessarily an infinite power series in its arguments, 
unlike any presently-known exact solution. However, it is apparently not trivial to show 
this, if in fact it is true. Complicating matters is that there are many possibilities for 
a base of size m, as m becomes large, and it is not clear that they will all predict the 
same polynomial, or reduce to the same function, K(p,r, s,t,u,v). Sorting out these 
issues would provide crucial insight into the question of whether or not thresholds of 
these problems are algebraic numbers. 

In this paper, I have found the approximate thresholds of the (4, 6, 12) and (3 4 , 6) 
lattices, completing the program begun in [22]. I also showed how this approximation 
leads to polynomial sequences for unsolved problems, and calculated the second 
polynomials for the (3 3 ,4 2 ) and (4,8 2 ) lattices. Finally, I showed how these results give 
some clues about whether or not generic percolation thresholds are algebraic numbers. 
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a) b) 

Figure 9. a) a base for the (3 3 ,4 2 ) lattice extending over m unit cells; b) the lattice 
resulting from setting pi — T{ = ti = for i > 1. The long red paths between the 
triangles can be contracted by setting Si = ti = 1 for i > 1 to recover the process 
whose base is a single unit cell. 
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Appendix 

Here I lay out the steps to finding the approximate thresholds of the (4, 6, 12) and 
(3 4 , 6) lattices. The procedures to find the other results in the paper are similar, and 
all thresholds are included in the supplemental material. 

(4, 6, 12) lattice 

First, we note the final threshold should have the form: 

FST(pi, p 2 ,p 3 ,p 4 ,p 5 , n,r 2 ,r 3 , r 4 , r 5 , si, s 2 , s 3 , s 4) s 5 ) = 

(1 - pi)L {p 3 ,p i ,p 5 ,r 1 ,r 2 ,r 3 ,r4,r 5 ,r 6 ,s 1 ,S2,t2,u 1 ,x,y) + 

PiLi(p 2 , p 3 , Pi, Ps, Pe, n,r 2 , r 3 , r 4 , r 5 , r 6 , Si, s 2 , h,t 2 , Wi, u 2 ) (20) 

with x = t\p 2 and y = u 2 p e . Starting with the L branch, all the lattices of which are 
in Figure [10] with the assignments of probabilities shown in Figure [121 we write 

L (P3, P4, P5, n, r 2 , r 3 , r 4 , r 5 , r 6 , s u s 2 , t 2 , u u x, y) = 
(1 - p 4 )L 00 (ri, r 2 , r 3 , r 4 , r 5 , r 6 , t 2 , Hi, u>, a;, y, z) + 

P4^oi(P3, Ps, r 2 , r 3 , r 4 , r 5 , r 6 , Si, s 2 , t 2 , Wi, x, y) (21) 

with w = s 2 p 3 and z = P5S1. Continuing, 

Looin, r 2 , r 3 , r 4 , r 5 , r 6 , £ 2 , u l5 x, y, z) = 
(1 - y)L 000 (a, n, r 4 , r 5 , r 6 , t 2 , r 2 r 3 , ty, x) + 

yLooi(ri, r 2 ,r 5 ,r 6 ,t 2 , 1 - (1 - r 3 )(l - ui),w,x,z). (22) 

L oo and L oi are actually the same lattice (Figure fTUlfd)) but with different probability 
labels. Thus, we have 

Looi{ri,r 2 ,r 5 ,r 6 ,t 2 ,v,w,x,z) = L 00 o(z 1 r ll v i r 5 ,rB 1 t2 1 r 2 ,w 1 x). (23) 

L oo can be written solely in terms of the result for the (4, 8 2 ) lattice in equation (jTJ), 

L 00 o(a, n, r 4 , r 5 , r 6 , t 2 , u, u>, x) = (1 - r 5 )FE(r 1( u, t 2 r 4 , xr 6 , w, a) + 
r 5 FE(r 1; t;,r 4 ,r 6 , W [l - (1 - x)(l - t 2 )], a) (24) 

and we have therefore determined L 00 . Turning to L i, we find that the A— lattice with 
some complicated bonds appears: 

L i(p 3 , Ps, ri, r 2 , r 3 , r 4 , r 5 , r 6 , si, s 2 , t 2 , m, x, y) = 
(1 - r 4 )A(gi,g 2 ,s 2 ,Si,ri) + 

uL 011 (pa, p 5 , n, r 2 , r 3 , r 5 , r 6 , si, s 2 , t 2 , Wi, x, y) (25) 
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with q x = p 3 r 6 [l - (1 - x)(l - r 5 t 2 )} and q 2 = p 5 r 2 [l - (1 - - r 3 «i)]. Next, 

Lqii(P3,Ps, r i, r 2 , r 3 , r 5 , r 6 , s x , s 2 , t 2 , «i, x, y) = 
(1 - r 3 )LQiio(p3,P5,r 1 ,r s ,r B ,s 1 ,82,t2,ux,yr 2 ,x) + 

r 3 L m (p 3 , n,r 2 , r s , r 6 , s 1; s 2 , t 2 , v, x) (26) 

with t> = p 5 [l — (1 — «i)(l — y)]. i^oiio can be expressed in terms of a lattice used in 
[22] and called the "rocket" lattice (Figure fTOTh); the full threshold was included in the 
supplementary material of [22] )• We have 

A)iio(P3, P5, r 5 , r 6 , si, s 2 , t 2 , Mi, x) = 
(1 - r 5 )R(si, s 2 , r u xr 6 , v, t 2 u 1 ,p 3 , p B ) + 

r 5 R(si,s 2 ,r 1 ,r 6 ,w,M 1 ,p 3 [l - (1 - a;)(l - t 2 )],p 5 ). (27) 

Now, 

A)in(P3, ri, r 2 , r 5 , r 6 , si, s 2 , i 2 , v, x) = 
(1 -r 5 )R(si,s 2 ,ri,a;r 6 ,r 2 ,t 2 ,p3,t;) + 

r5-^onu(si, s 2 ,r 2 ,r e ,ri,p 3 [l - (1 - x)(l - £ 2 )],t>) (28) 
and Lomi can either be expressed in terms of the rocket lattice or 

L uii{p,r,s,t,u,v,x) = (1 - w)H([l - (1 - pt){\ - rs)],v,x) + 

wH([l - (1 - p)(l - r)] [1 - (1 - t)(l - s)],v, x). (29) 

This completes the Lq branch, so now we need to find L\. The lattices needed for this 
are shown in Figure [TT] with probability assignments in Figure fT3l To begin, 

£i(P2,P3, P4, P5, Pe, r x , r 2 , r 3 , r 4 , r 5 , r 6 , s l} s 2 , * 1? t 2 , «i, w 2 ) = 
r 4 ^n (p 2 ,P3, P4, Ps, Pe, rl r 2 , r 3 , r 5 , r 6 , s x , s 2 , t 1} t 2 , u u u 2 ) + 
(1 - r 4 )L 10 (p 2 ,p 3 ,p A ,p 5 ,p 6 ,r 1 ,r 2 ,r 6 ,s 1 ,s 2 ,t 1 ,u 2 ,r 3 u 1 ,r 5 t 2 ), (30) 

^10(P2,P3,P4,P5,P6,n,^2,^6,Sl,S 2 ,tl,M 2 ,f,w) = 

riLoii(r 6 ,r 2 ,p 4 ,p 5 ,p 6 ,p 2 ,p 3 ,s 2 ,Si,ti,M 2 ,w,w) + 

(1 - ri)Lioo(p2,P3,P4,P5,P6,^i,M2,w,w,Sir 6 ,s 2 r 2 ) (31) 

and 

Lloo(p2,P3,P4,P5,P6,h,U 2 ,V,W,X,y) = 

(1 - t 1 )L 100 o(P2,P3,P4,P5,P6,U 2 ,V,W,X,y) + 

£i£iooi(P4,P5,P6,W2,^,a:,2/,P3[l - (1 - tw)(l -P2)])- (32) 
Liooi can be found by 

Liooi(P4,P5,P6,u 2 ,v,x,y,z) = (1 - v)A(x,yu 2 ,z,p 5 p 6 ,p 4 ) + 

vFE(l - (1 -p 6 )(l - u 2 ),z,pi,p 5 ,x,y). (33) 
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L 1000 is similarly determined: 

(1 - u 2 )FE(p 6 p 2 ,p 3 ,p 4 ,p 5 ,z,yv) + 
u 2 FE(p 2 ,p 3 ,p A ,p 5 [l - (1 -v)(l -p 6 )],z,y). 

Moving on to Lu we find 

Ln(p 2 , p 3 , Pa, Ps, Pe, n, r 2 , r 3 , r 5 , r 6 , Si, s 2 , h,t 2 , m, u 2 ) = 
(1 - r 2 )L 11 o{p 2 ,P3,P4,p 5 ,P6,r 5 ,r 6 ,si,t 1 ,t 2 ,u 1 ,r 1 s 2 ,r 3 u 2 ) + 
r 2 L lu (p 2 ,p 3 ,p 4 ,p 5 ,p (i ,r 1 ,r 3 ,r 5 ,r 6 ,s 1 ,s 2 ,t 1 ,t 2 ,u 1 ,u 2 ), 

and 
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(34) 



(35) 



Lni(p 2 ,P3, P4, P5, Pe, n, r 3 , r 5 , r 6 , s x , s 2 , ti, t 2 , u u u 2 ) = 

(1 - r 6 )L ino (p 2 ,p 3 ,p 4 ,p 5 ,p 6 ,r 3 , s 2 ,t 2 ,m,u 2 , sir^ttrs) + 

r 6 L m i (p 2 , p 3 , pi, p 5 , pe,n, r 3 , r 5 , s x , s 2 , ti,t 2 ,U!, u 2 ). (36) 

The lattice L 110 can be obtained by the threshold for L that we have already derived: 

Luo(P2,P3,P4,P5,P6,r5,re,s 1 ,tx,t 2 ,u 1 ,v,w) = 

and similarly, 

^lllo(P2,P3,P4,P5,P6,r3,S2,t 2 ,'Ul,W2,f,w) = 
L (l,r 3 , l,p 6 , l,p 2 ,P3,P4,P5,Ml,M2,S 2 ,t 2 ,W,w). 



(37) 



(38) 



Further, 



where 



Lun(p2,P3,PA,P5,P6,n,r 3 ,r 5 , si, s 2 ,*i,*2, «i, w 2 ) = 
(1 - p 5 )^nno(P2,P3, ri, r 3 , r 5 , s 2 , t ls t 2 , u 2 , itip 6 , sip 4 ) + 
Ps-^iiin(P2, P3, P4, Pe, n, r 3 , r 5 , si, s 2 , h,t 2 , u x , u 2 ) 

Luuo(p2,P3, n,r 3 , r 5 , s 2 , h,t 2 , u 2 , v, w) = 
L (l,p 2 ,p 3 ,r 5 , l,n, l,r 3 , l,t 2 ,ti,M 2 ,s 2 ,f,w) 



(39) 



(40) 



and 



^11111 (P2, P3, P4, P6, r l, r 3, 5"5, S l> S 2, *1, *2, «1, « 2 ) = 

(1 - p 3 )L imi0 (p 6 , n, r 3 , r 5 , si, ti, Mi, U2,p 2 t 2 ,p 4 s 2 ) + 

P3Luiiu{P2,Pi,P6, n,r 3 , r 5 , s x , s 2 , t 1} t 2 , u t , u 2 ). (41) 

The £mm lattice is just the Archimedean (3,4,6,4) lattice (Figure [TT^l)), for which 
the linear threshold, denoted TFSF, was found in [22]. Therefore, 

£111111 (P2,P4, Pe, n,r 3 , r 5 , Si, s 2 , ti, t 2 , «i, w 2 ) = 

TFSF(si, s 2 , n,p 4 , r 3 , w 2 , «i,p 6 , ti, r 5 , t 2 ,p 2 ). (42) 



Polynomial sequences for bond percolation critical thresholds 



18 





g) 



h) 







Figure 10. The lattices from the "O-branch" in the (4,6,12) tree, a) Lq\ b) Loo! 
c)L i; d) L om and L 001 ; e) Lou 5 f) ^ono! g) loxil! h) rocket lattice; g) Iqiiii- 



£111110 can be found from a previous result, 

£iimo(P6, ri, r 3 , r 5 , si, £i, «i, u 2 , v, iw) = 

•^mio (pe, 1, n>, n, r 3 , ti, Mi, w 2 , si, w, u). (43) 

We have now determined all the sublattices, and we can climb back up the tree to find 
the linear threshold for the (4, 6, 12) lattice. The full expression has 1932 terms and is 
included in a text file in the supplementary material. 

(3 4 , 6) lattice 

The lattices in this tree are shown in Figures HH and [151 except the (3, 4, 6, 4) lattice 
(TFSF), which is in Figure [TTT l) . Probability assignments are in Figures [T6l and [T71 
Denoting the threshold TFS, we have 

TFS(pi, p 2 , P3,P4, Pb, n, r 2 , r 3 , r 4 , r 5 , s u s 2 , s 3 , s 4 , s 5 ) = 

(1 - p 5 )L (pi,p 2 , Pa, P4, n, r 2 , r 3 , r 4 , r 5 , s u s 2 , s 3 , s 4 , s 5 ) + 

p 5 Li(?T, r 2 , r 3 , r 4 , r 5 , Si, s 2 , s 3 , s 4 , s 5 , t, w) (44) 
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Figure 11. The lattices from the "1-branch" in the (4,6,12) tree, a) L\\ b) L w ; 
c ) Lu; d) iioo; e) Lu ; f) Lin; g) Liooo; h) Liooi; i) Lino; j) Lim; k) Liimo; 1) 
(3,4,6,4) lattice. 



with t — 1 — (1 — Pi)(l — P3) and w = 1 — (1 — p 2 )(l — ^4)- Tackling the 0-branch first, 
we have 

Lo(pi,P2, P3, Pa, n, r 2 , r 3 , r 4 , r 5 , si, s 2 , £3, s 4 , s 5 ) = 
(1 - S5)-^oo(pi,P2,P3,P4,ri,r2,r3,r4,r 5 ,si,S2, 53,54) + 
s 5 ^oi (pi, P2, P3, P4, ?T, r 2 , r 3 , r 4 , r 5 , it, u) (45) 
where u — 1 — (1 — Si)(l — S4) and t> = 1 — (1 — s 2 )(l — S3), 

£01 (pi, P2, Pa, P4, ?T, r 2 , r 3 , r 4 , r 5 , u, u) = 

(1 - u)Loio(pi,P2,P3,P4,ri,r2,r 3 ,r4,r 5 ,u) + 

u£on(Pi,P2,P3, ?T, r 2 , r 3 , r 5 , w), (46) 
with w = 1 — (1 — r 4 )(l — r 5 ), 

^oii(Pi,P2,P3,^i,^2,r 3 ,r 5 ,v,w) = 

PiL 0111 (p 1 ,P3 : r 1 ,r2,r 3 ,r 5: v,w) + (1 - Pi)£ iio(Pi,P3, ?T, r 3 , r 5 , u, iu), (47) 
^ono(pi,P3,n,r2,r3,r 5 ,u,«;) = 
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Figure 12. Unit cells with probability assignments for lattices from the "0-branch" in 
the (4,6, 12) tree, a) L ; b) L oa ; c)Lq\\ d) i oo and L Q01 ; e) £ n ; f) Anio; g) Ami! 
h) rocket lattice; i) Lqiiii- 



(1 - r 3 )DS(wpi,r 2 ,r 5 ,p3,r 1 ) + r 3 L iioi(pi,P3, »" 2 , v, x), (48) 

where x = 1 — (1 — w)(l — r^) and DS denotes the "decorated square" lattice shown in 
Figure IT4T n) with the threshold that was given in [22] , 

DS(p, r, s, t,u,v) — 1 — p — st — ru — rtv — suv + rstv + 

rstu + rsuv + rtuv + stuv — 2rstuv. (49) 

Continuing, 

Amoi(Pi,P3,?"i,r2,'u,x) = (1 - r 1 )B(p 1 ,p 3 r 2 ,x,v) + 
nHfa, 1 - (1 - x)(l - r 2 ), 1 - (1 - v)(l - pa)), (50) 
where B denotes the martini-B lattice, 

L ni(pi,P3,ri,r 2 ,r 3 ,r 5 ,v, w) = (1 - w)FED(w, n, r 5 , r 3 pi,p 3 , r 2 ) + 
FED(w, n, 1 - (1 - r 5 )(l - r 3 ), 0, 1 - (1 - p 4 )(l - px), r 2 ), (51) 
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Figure 13. Unit cells with probability assignments for lattices from the "1-branch" 
in the (4, 6, 12) tree, a) L\\ b) L w ; c) L n ; d) L wo ; e) L no ; f) imi g) £1000; h) Liooi; 
i) Lino; j) lam; k) Liimo; 1) (3,4,6,4) lattice. 



where FED is the dual of the (4, 8 2 ) lattice (Figure fl47o)). The critical surface for a 
dual graph is generically derived from the critical surface of the original graph by the 
substitutions Pi — > 1 — p^, and multiplication by a minus sign to ensure the constant 
term is equal to 1. Thus, in the (4, 8 2 ) case we have, 

FEDQo, r, s, t, u, v) = -FE(1 - p, 1 - r, 1 - s, 1 - t, 1 - u, 1 - v). (52) 

Now, 

Loo(pi,P2,P3,Pa, ru r 2 , r 3 , r 4 , r 5 , s ls s 2 , s 3 , s 4 ) = 

r 5 Looi(pi,P2,P3,P4,Si,S2,S3,S 4 ,M,w) + 

TFSF(pi,p 2 ,P3,P4, ri, r 2 , r 3 , r 4 , s x , s 2 , s 3 , s 4 ), (53) 
with it = 1 — (1 — J"i)(l — r 2 ), t> = 1 — (1 — r 3 )(l — r 4 ), 
^001 (Pl,P2,P3, Pi, si, s 2 , s 3 , s 4 , it, u) = 
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(1 - p,i)L Q010 (p 2 ,p 3 , s 1 , s 2 , s 3 , s 4 ,u, x) + 

P4L 0011 (p 1 ,p 2 ,P3,S 1 ,S2,S 3 ,U,w), (54) 

with x = vp\ and w — 1 — (1 — t>)(l — s 4 ), 

^001o(P2,P3,Sl,S2,S 3 ,S 4 ,M,x) = 

P3-^0010l(P2, Si, S3, S A ,X,y) + (1 -p 3 )KD(Mp 2 ,S3,S4,X, Si,S 2 ), (55) 

where y = 1 — (1 — u)(l — s 2 ) and KD denotes the kagome dual or dice lattice, shown 
in Figure fT4Tp). 

L in(pi,P3,r 1 ,r 2 ,r 3 ,r 5 ,v,w) = (1 - w)FED(v,r 1 ,r 5 ,r 3 p 1 ,p 3 ,r 2 ) + 
wFED{v,r u 1 - (1 -r 5 )(l -r 3 ),0,l - (1 - p 3 )(l - p x ), r 2 ), (56) 

^010(^1,^2,^3,^4, ri ) r2,r 3 ,r4,r 5 ,w) = 
P4-^oioi(pi,P2,P3, n, r 2 , r 3 , r 4 , r 5 ) + 

(1 -p4)£oioo(pi,P2,P3,n,r 2 ,r 3 ,r 4 ,r 5 ,t;), (57) 
^oioo(Pi,P2,P3, ri, r 2 , r 3 , r 4 , r 5 , t>) = 

r 5 A(pi[l - (1 - r)(l - r 4 )],p 2 ,p 3 ,v, 1 - (1 - ri)(l -r 2 )) + 

(1 - r 5 )L 010 oo(pi,p 2 ,p 3 ,r 1 ,r 3 ,M,f), (58) 

where u = r 2 r 4 , 

L iooo{pi,P2,P3,n,r 3 ,u,v) = 

nH(p 2 , 1 - (1 - p 3 )(l - v), Pl [l - (1 - r 3 )(l - it)] ) + 

(1 - ri)A(pi,p 2 ,p 3 ,ur 3 ,w), (59) 

^oioi(Pi,P2,P3, n, r 2 , r 3) r 4 , r 5 , v ) = 

(1 - r 5 )L iio(l, r 3 , n, u, p 2 , p 3 , r 2 r 4 , p x ) + 

r 5 L no(l, 1 - (1 - r 4 )(l - r 3 ), 1 - (1 - n)(l - r 2 ), v,p 2 ,p 3 , 0,pi), (60) 
^001l(Pl,P2,P3,Sl,S 2 ,S 3 ,M,w) = L io(s 3 ,si,s 2 ,w,p 3 ,pi,p 2 , 1,0, u), (61) 

and finally, 

£00101 (P2, si, s 3 , s 4 ,x,y) = L 010 (s 1 ,s 3 ,s 4 ,y,0,p 2 ,l,l,0,x) (62) 
This takes care of the 0-branch. Moving on to the 1-branch, we have 

Li(ri, r 2 , r 3 , r 4 , r 5 , s x , s 2 , s 3 , s 4 , s 5 , t, w) = 

(1 - t)Li (ri, r 2 , r 3 , r 4 , r 5 , si, s 2 , s 3 , s 4 , s 5 , w) + 

tL u (r 2 , r 3 , r 4 , r 5 , s x , s 3 , s 4 , s 5 , w, (63) 
where v = 1 — (1 — s 2 )(l — ri), 

Lu(r 2 , r 3 , r 4 , r 5 , si, s 3 , s 4 , s 5 , w, = 
r 2 £ui (r 3 , r 4 , Si, s 3 , s 4 , s 5 , w, w) + 

(1 - r 2 )Luo(r 3 , r 4 , r 5 , si, s 3 , s 4 , s 5 , u, v) (64) 
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with w = 1 — (1 — v)(l — r 5 ), 

£111 (r 3 , r 4 , si, s 3 , s 4 , s 5 , u, w) = 

(1 - u)Lnio(r 3 , r 4 , si, s 3 , s 4 , s 5 , u>) + 

uMD(si,t,s 5 ,iw,S3,r 3 ) (65) 

where t = 1 — (1 — r 4 )(l — s 4 ) and MD denotes the dual of the martini lattice (equation 
©). Now, 

£iiio0"3, ?*4, s h s 3 , s 4 , s 5 , w) = 

S(l - (1 - si)(l - s 4 s 5 ), 1 - (1 - r 4 )(l - w 3 )) + (66) 
s 3 AD(r 4 ,si,s 4 ,r 3 , 1 - (1 - to)(l - s 5 )) (67) 

where S represents the square lattice, and AD is the dual to the martini-A lattice. Now, 

£iio(t* 3 , r 4 , r 5 , si, s 3 , s 4 , s 5 , w, t>) = 

(1 - w)Lnoo(r 3 , r 4 , r 5 , si, s 3 , s 4 , s 5 , u) + 

^1101(^3,^4, r 5 , s 1 ,s 4 ,u,w) (68) 
with w = 1 — (1 — Ss)(l — s 3 ), 

£llOo( r 3, r 4, r 5? s l, s 3; s 4, S5, u) = 

(1 - s 5 )Lnooo(^3, n, r 5 , si, s 3 , s 4 , u) + 

55-^11001 (t* 3 , 7*4,7*5, S3, U,V) (69) 

where u = 1 — (1 — si)(l — s 4 ), 

£11001 fo, r 4 , r 5 , s 3 , u, u) = L mo (s 3 , 0, u, r 5 , v, r 4 , r 3 ), (70) 

and, 

£11000(73, 7*4, r 5 , si, s 3 , s 4 , u) = (1 - r 4 )FE(s 4 , w, r 3 , r 5 , s 3 , s x ) 

r 4 S(s 3 [l - (1 - w)(l - s 4 )][l - (1 - r 5 )(l - r 3 )], s x ). (71) 

Continuing, 

£ll0l(r3, 7"4, 7- 5 , Si, s 4 , m, w) = 

s 4 B(l - (1 - si)(l - w), r 3 , 1 - (1 - u)(l - r 4 ), r 5 ) + 

(1 - s 4 )AD(r 3 ,ww, si,r 5 ,r 4 ), (72) 

£io(?"i, r 2 , r 3 , r 4 , r 5 , si, s 2 , s 3 , s 4 , s 5 , u) = 

(1 - ri)Li 00 (r 2 , r 3 , r 4 , r 5 , s 4 , s 2 , s 3 , s 5 , s 5 , u) + 

ri^ioi (r 3 ,r 4 , si, s 2 , s 3 , s 4 , s 5 , w,v), (73) 

where t> = 1 — (1 — r 2 )(l — r 5 ). L 101 is actually the same lattice as L0101, but with 
different labels for probabilities. It can also be written in terms of £110, 

£101(7-3, 7*4, Si,S 2 , s 3 ,s 4 , s 5} u,v) = 

(1 - s 5 )Liio(r 3 , r 4 , v, s x s 2 , s 3 , s 4 , 0, u, 1) + 

S5£no(7* 3) r 4 ,f,0, s 3 , 1 - (1 - si)(l - s 4 ), s 2 ,u, 1). (74) 
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Now, 

-^100(^2,^3,^4,^5, Si, s 2 , s 3 , s 4 , s 5 ,w) = 

(1 - r 5 )L 01 oo(r 3 , r 2 , r 4 , s 4 , si, s 3 , s 2 , s 5 , w) + 

^5^1001(^2, si,s 2 , s 3 , s 4 , s 5 ,w, 1 - (1 -r 3 )(l -r 4 )), 

^1001(^2, Si, s 2 , S3, s 4 , s 5 , it, v) = 
s 5 A(r 2 , 1 - (1 - s 2 )(l - s 3 ), 1 - (1 - si)(l - s 4 ), u, 
(1 - S 5 )Liooio(^2, s 3 , s 4 , it, iw) 
with w = sxs 2 . At long last, 

£10010 fa, s 3 , s 4 , it, iu) = (1 - s 4 )B(r 2 , w, v, s 3 u) + 
s 4 S(r 2 [l - (1 - u)(l - v)\, 1 - (1 - w)(l - s 3 )). 
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Figure 14. The lattices from the "0-branch" in the (3 4 ,6) tree, a) L ; b) L i; c)Loo; 
d) ^ooi; e) Loio ; f) ion; g) ^oioo; h) Loioi! i) -^ono; j) Loin; k) Looio; 1) Loon; m ) 
£00101; n ) decorated square; o) (4,8 2 ) dual; p) dice lattice (kagome dual) 
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Figure 15. The lattices from the "1-branch" in the (3 4 ,6) tree, a) L\\ b) L w ; c)In; 
d) Lioo; e) Lno ; f) £111; g) £1110; h) ^1001; i) iiioo! j) ^lioi; k) Lnooi; 1) £11000 
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Figure 16. Unit cells with probability assignments for lattices from the "0-branch" 
in the (3 4 ,6) tree, a) L ; b) L i; c)L Q0 ; d) L a0 i; e) £ io ; f) L on ; g) £ ioo; h) £0101; 
i) L 01W ; j) Lom; k) iooio; 1) ^oou; m) Looioi; n ) (4, 8 2 ) dual; o) dice lattice (kagome 
dual) 




Figure 17. Unit cells with probability assignments for lattices from the "1-branch" 
in the (3 4 ,6) tree, a) L\\ b) L w ; c)L n ; d) L 100 ; e) L no ; f) L m ; g) i mo ; h) iiooi; 
i) iiioo; j) inoi; k) Lnooi; 1) £11000 



